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ABSTRACT 



Context. The study of the dust extinction in high-redshift galaxies is fundamental to obtain an estimate of the corrected Star Formation 
Rate Density (SFRD) and to put constraints on galaxy evolution models. 

Aims. We plan to analyse dust extinction in Lyman Break Galaxies (LBGs) by introducing a new and more reliable approach to 
their selection and to the characterization of their distribution of UV slopes /?, using deep IR images from HST We fully discuss the 
methodology and the results focusing on a robust sample of z~4 LBGs. 

Methods. We exploit deep WFC3 IR observations of the ERS and HUDF fields over GOODS-South, combined with HST-ACS 
optical data, to select z~4 LBGs through a new (B-V) vs. (V-H) colour diagram. The UV slope of the selected galaxies is robustly 
determined by a linear fit over their observed I, Z, Y, J magnitudes, coherently with the original definition of /?. The same fit is used to 
determine their rest-frame UV magnitudes M l600 through a simple interpolation. We estimate the effect of observational uncertainties 
with detailed simulations that we also exploit, under a parametric maximum-likelihood approach, to constrain the probability density 
function of UV slopes PDF(/3) as a function of rest-frame magnitude. 

Results. We find 142 and 25 robust LBGs in the ERS and HUDF fields respectively, limiting our sample to S/N(H)>10 objects. Our 
newly defined criteria improve the selection of z ~4 LBGs and allow us to exclude red interlopers at lower redshift, especially z~3-3.5 
objects. We show that the use of a linear fit to estimate ft and an accurate characterization of observational effects are required in 
such analysis of flux limited samples. We find that z ~4 LBGs are characterized by blue UV slopes, suggesting a low dust extinction: 
all L < L* galaxies have an average UV slope (J3) =s —2.1, while brighter objects only are slightly redder {{f}) m -1.9). We find an 
intrinsic dispersion ^ 0.3 for PDF(/3) at all magnitudes. The SFRD at z~4 corrected according to these estimates turns out to be lower 
than previously found: log(SFRD)=: -1.09 M Q /yr/Mpc 3 . Finally, we discuss how the UV slope of z~4 galaxies changes as a function 
of the dust-corrected UV magnitude (i.e. SFR). We show that most galaxies with a high SFR (> 80 M Q /yr) are highly extincted 
objects. Among galaxies with lower SFR, we detect many with a much lower amount of reddening, although current observational 
limits prevent us from detecting those with high extinction, if they exist. 

Key words. Galaxies:distances and redshift - Galaxies: evolution - Galaxies: high redshift 
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1. Introduction 

It is now nearly 20 years since the first surveys explic- 
itly de dicated to detect galaxies at z — 3 have been com- 
pleted (ISteidel & Hamilton! 1 19921 [l993t iGiavalisco et~allll994t 
Stei del et al.ll 19951 ). In the early 90's, these seminal works have 
introduced the so-called "Lyman Break" method to select rest- 
frame UV bright objects, i.e. star-forming, moderately absorbed 
galaxies, at z - 3 and at higher redshifts. In the following 
years, this technique has been widely appli ed by surveys cov- 
erin g either progressively wider areas (e.g. LSteidel e"tai1l2003t 
lOuchi et all l2004t iMcLure et all 120091) to enla rge the volume 
samp l ed, or reaching fainter magnitudes (e.g. Bouwens et alj 
l2007t iReddv & Steide]| l2009) to detect the faintest galaxies, or 
extending to longer wavelengths to explore the highest redshifts 
(e.g.lBouwens et al.l2004tlCastellano et alj2010HBouwens et alj 
1201 11) . In the meantime, spectroscopic follow-up campaigns 
of Lyman Break Galaxies (LBGs) samples have shown that 
this technique is extremely efficient and reliable up to the 



highest redshifts (e.g. ISteidel et al.lll996t [Vanzella et alj [2 009: 
Star k et alJl20Tot IVanzella et alJl201 lHPentericci et al.ll201 ll) . 



Some of the primary goals of these surveys has been the sta- 
tistical analysis of high redshift galaxies, in particular of their 
UV Luminosity Fu n ction (LF) from z —3 to z -7 and beyond 
(ISteidel et aljl 19991: iBouwens et all 120071: IMcLure et all [2009; 



iBunker et al.ll20l6T iGrazian et al. 20111) . and their contribution 
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to the Star Formation Rate Density (SFRD). 

Understanding the physical properties of UV-selected galax- 
ies is also fundamental to derive the global properties of the 
Universe at these high redshifts. In particular, to estimate their 
contribution to the total SFRD one needs to apply the proper con- 
version between the observed UV rest frame luminosity and the 
ongoing Star Formation Rate (SFR). This conversion factor de- 
pends on the physical properties of the stellar population (IMF, 
metallicities and ages) and on the amount of dust extinction. The 
slope of the UV continuum in LBGs provides precious informa- 
tion to measure these quantities. Apart from weak absorption 
lines and other minor features, the spectrum of a star-forming 
galaxy longward of the Lyo- wavelength can be reasonably ap- 
proximated by a power law until 3200A for most combinations 
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of stellar properties and attenuation curv es dCaIzettietalJ[ l994: 
lMeureret al.1 119991: ICalzetti et al.l l2000b iGallerani et alT feOlO). 
The actual slope is the result of several factors involved in the 
stellar mixture in the galaxy: the total amount and composi- 
tion of dust grains, the metallicity Z, the initial mass function 
(IMF), and the star formation history (SFH). Lower metallic- 
ity stars, top heavy IMFs and young ages all tend to produce 
bluer spectra than higher metallicities, standard IMFs and aging 
stellar populations. Nebular continuum emission can also con- 
tribute to redden the slope for galaxies. We note that nebular 
emission is expected only if the escape fraction of UV ionizing 
photons is small, as currently observed in most of t he galaxies 
at z < 3.5 (IVanzella et al.ll2010fc iBoutsia et al.ll201 ll) . However, 
dust absorption is likely the strongest reddening factor, and the 
observed slope is often directly converted into dust extinction 
assuming standard stellar populations. 

In all cases, the UV spectrum is well approximated by a sim- 
ple power law fx = AP: the accurate determination of the ex- 
ponent /3 is the central goal of the present paper. For reference, 
a dust-free stellar population with solar metallicity and constant 
star-formation rate has /3 - -2.2. Clearly, it is impossible to con- 
strain all the physical quantities mentioned above from a single 
parameter. Nevertheless, a robust measure of (3 is the first manda- 
tory step to constrain these free parameters and derive the global 
properties of LBGs. 

Despite the large samples collected in the last years, only a 
few estimates of the statistical distribution of f3 are available to- 
day, and no firm conclusion has been reached on its evolution 
with luminosit y and redshift. Previ ous determinations have been 
provided bv lBouwens et alJ ([2009), that found a strong decrease 
of the average UV slope both wit h rest-frame magnitude and 
redshift from z ~ 2 to z ~ 6, and bv lBouwens et alJ (|2010c) pre- 
senting e vidence for very blue UV slopes in their z~7 sample. 
At z~6-8 Finkelstein et al. (2010) found blue UV slopes con- 
sisten t with dust-fr e e nor mal stellar populations. On the other 
hand, iDunlop et al.l d201 ll) found that UV slopes are consistent 
with /? =; -2 at z ~ 5 - 7, and argue that very steep values found 
at high redshif t are due to the effe cts of photometric scatter and 
uncertainties. IWilkins et al] (120111) reached a similar conclusion 
from their z ~7 sample, but at z~ 5 fo und average slopes red - 
der than both | Dunlop et al] d201ll) and iBouwens et al] d2009). 
At z ~ 4 lOuchi et all (l2004l) found an average E(B-V)^ 0.15 and 
no si gnificant depende nce on UV magnitude for L > L» LBGs, 
while Lee et a i] (l20Tl found a decreasing extinction, from E(B- 
V>0.22 at L =s 4 - 5L» to E(B-V)^0.07 at L * L„. 

Discrepancies among different works are prob ably due to the 
know n limits of LBG selection techniques (e.g. IStanwav et alJ 
2008) and to uncertainties in the measure of the UV slope, of- 
ten based on a single colour. These limitations can be overcome 
by assembling deep near IR observations of large LBG samples: 
these are required to strengthen colour selections and to fully 
sample the UV continuum up to 3000A providing a proper 
estimate of /3. For this reason, early results based on optically 
selected samples are intrinsically prone to larger systematic ef- 
fects. Even more importantly, however, they suffer of a concep- 
tual limitation. For a given dereddened luminosity (or equiva- 
lently a given SFR), even a modest dust absorption produce a 
significant dimming of the galaxy, which may escape detection 
in optically selected surveys. 

In this paper we take advantage of the new IR observations 
carried out with WFC3 on board of HST. We show that these 
deep observations in the Y, J and H bands provide significant 
advantages over previous ones in three different aspects: 



• They make possible a much cleaner selection of LBGs, 
compared to the traditional criteria that use optical bands only. 
This allows us to reject lower redshift interlopers with much 
higher reliability. 

• They yield a robust determination of the UV slope, provid- 
ing both a large leverage in wavelength and a sizeable number 
of independent determinations of the flux (of each object) at dif- 
ferent wavelengths. 

• They allow us to detect galaxies with a larger slope for 
a given dereddened luminosity. In this respect, the H-band that 
we use to detect galaxies constitutes the optimal choice since 
it corresponds to the reddest portion of the UV rest-frame do- 
main (3200A rest-frame at z - 4), significantly redder than the 
1500 - 1800A of optically selected samples. Hence our samples 
are intrinsically more complete in terms of content of reddened 
galaxies than optically selected ones, although this is currently 
limited by the current depth of optical observations used for the 
colour selection. 

To fully exploit these potential advantages, we apply a tai- 
lored colour selection criterion to our H-band selected cata- 
logue, and we measure /3 with a straightforward linear fit to the 
observed magnitudes. In principle, similar results could be de- 
rived usin g a SEP fitting appr oach on the same H-selected sam- 
ple (e.g. iFontana et al]|2003l) . where the slope depends on the 
free parameters mentioned above: IMF, metallicity, age, E(B-V). 
We decide to adopt a simpler method for two reasons. First, as 
mentioned above, the four free parameters are to a large extent 
degenerate, even using additional constraints from the Spitzer 
magnitudes. Our approach determines a single parameter (J3) 
from four observed magnitudes. Further physical interpretation 
can be derived from the observed values of f3 alone. In addition, 
a complete analysis requires extensive Monte Carlo simulations, 
that are required to take into account the involved systematic ef- 
fects. Selecting galaxies with a simple colour selection criterion 
and adopting a single free parameter makes these simulations 
affordable. 

In this paper we apply our method to z - 4 LBGs. In fu- 
ture papers, we plan to extend the analysis to higher redshifts, 
using a lso the forthcoming data sets provided by th e CANDELS 
survey dGrogin et alj|201 ltlKoekemoer et ai]|201 ll) . 

The plan of the paper is the following: in Sect. [2] we present 
our data-set from which we select a sample of B-dropout galax- 
ies applying a newly defined IR-based colour selection, dis- 
cussed in Sect. [3] In Sect. [4] we present the estimate of UV 
slopes and rest-frame magnitudes for our objects, along with the 
simulations used to evaluate systematic effects and uncertain- 
ties. These simulations are exploited in Sect. [5] to constrain the 
distribution of /3 at different UV magnitudes, from which a dust- 
corrected value of the SFRD is estimated in a straightforward 
way (Sect. [6]). In Sect. [7] we discuss how the UV slope of z~4 
galaxies in our sample changes as a function of their dereddened 
UV emission. Finally, a summary is given in Sect. [8] 

Th roughout the whole paper, we adopt a ICalzetti et alJ 
(2000) extinction law, unless otherwise specified. Observed and 
rest-frame magnitudes are in the AB system, and we adopt the 
A-CDM concordance model (Hq = 70km/ s/Mpc, Q. M = 0.3, 
and Q A = 0.7). 



2. Data and Photometry 

In the present pape r we exploit data ob tained in the WFC3/IR 
ERS2 program dWindhorst et all 1201 lb and the ACS and 
WFC3/IR UDF programs dBeckwith et all2006HBouwens et alJ 
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Fig. 1. The (B-V) vs (V-H) colour selection window adopted in the present paper (right) along with the commonly adopted one 
based on the (V-Z) colour (left). Blue filled squares and lines are for a CB07 LBG model with Z/Z = 0.2 and Age=100Myr at 
increasing redshift and E(B-V) respectively. Models having solar metallicity and Age=100Myr and Age=300Myr are represented 
in magenta (crosses and long-dashed lines), and red (empty squares and short-dashed lines), respectively. Red and green points 
indicate lower redshift passively evolving and dusty star-forming galaxies respectively. See main body for details. 



I2010ah . w here the data were proc essed and combined as de- 
scribed in iKoekemoer et al.l d201 ll) . along with the full public 
ACS data available in the GOODS-South field (Giavalisc o et al.1 



2004). Full details on image red uctio n and catalogue bu ilding 
are given in iGrazian et al.l ( 1201 ll) and ISantini et al.l (1201 ll) : we 
summarize here the main properties of our dataset. 

The GOODS-ERS WFC3/IR dataset (HST Program ID 
1 1359) covers ~ AOsq.arcmin. in the Fo98> ^125, #160 filters. It is 
made of 10 contiguous pointings in the GOODS-South field ob- 
served during 60 HST orbits, for a total of 4800-5400s per point- 
ing and filter. The astrometry of WFC3 ERS images w as bound 
to the GOODS-S v2.0 z band dGiavalisco et al.ll2004l) . The lim- 
iting magnitudes are Y=27.3, J=27.4, and H=27.4 at 5<x in an 
area of ~ 0.1 1 arcsec 1 (corresponding to 2FWHM). 

The HUDF WFC3/IR dataset (HST Program ID 1 1563) is a 
total of 60 HST orbits in a single pointi ng of 4.7 sq.arcmin. (e.g. 
lOesch et al.l2010tlBouwens et al.l2010bl) in three broad-band fil- 
ters (16 orbits in Y\os, 16 in /125, and 28 in //i6o)- The HUDF 
images were reduced in the same way as the ERS ones and 
they have been astrometrically matched using the ACS HUDF 
sources. For the Y\q$ filter, we have selected a subset of 7 vis- 
its out of 9 that were not badly affected by persistence from 
bright sources. The magnitude limits of HUDF WFC3 observa- 
tions are Y=29.0, J=29.2, and H=29.2 total magnitudes for point 
like sources at 5<x in an aperture of ~ 0.1 1 arcsec 1 . 

In both fields we use the available BViz ACS public data: 
the latest V2.0 version of ACS BViz images of GOODS-South 
released by the STSci (M. Giavalisco and the GOODS Team, 
in prepar ation) for the ERS, and the deep BViz ACS data pre- 
sented in Bec kwith et al.l (120061) for HUDF. To ensure precise 
colours, the ACS bands were smoothed with an appropriate ker- 
nel to match the resolution of the WFC3 images. 

We perform object detection on the H160 band, obtaining 
10946 and 2418 sources for the ERS and HUDF, respectively. 
Their total magnitudes have been c omputed using MAG_BEST of 
SExtractor dBertin & Arn outs 1996) for galaxies more extended 
than 0.1 1 arcsec 1 , while circular aperture photometry (diameter 
equal to 2 times the WFC3 FWHM) is used for smaller sources. 
The magnitudes which are computed in a circular aperture have 
been corrected to total by applying an aperture correction of 0.4 



magnitudes. Colours in BVIZYJ have been measured running 
SExtractor in dual image mode, using isophotal magnitudes for 
all the galaxies. 

3. A new approach to the selection of z~4 LBGs 

We define the sample of z~4 galaxies by applying the Lyman- 
break technique. The main spectral feature that enables the iden- 
tification of high-z galaxies is the drop shortward of the Lyman- 
a, which, at z>3.5, is redshifted at A > 5500A: the Lyman break 
can thus be effectively sampled by the B and V bands (B-dropout 
galaxies). The usual colour-selection of z~4 LBGs is carried out 
on a (B-V) versus (V-Z) diagram (e g. Giavalisco et al.l [2004; 
iBouwenset all 120091: IVanzella et at] 120091) to separate LBGs 
from lower redshift Balmer break and dusty galaxies which show 
re dder colour at longer w avelengths. However, as already noted 
by iBeckwith et al.1 (120061) the redshift selection range provided 
by these criteria has a non-neglig ible dependence on the amount 
of extinction. For example. IBeckwith et al.l (12006) (Tab. 8 and 9) 
show that a lOOMyr star-forming model with E(B-V)=0.15 can 
be included in the selection from lower redshifts (z„„„ = 3.3) 
than a dust-free one (z m ,„ = 3.5). 



3.1. Colour selection 

In this paper, thanks to the very deep WFC3 near infrared im- 
ages, we adopt a modified version of the B-dropout selection 
which is based on the (V-H) colour instead of the ( V-Z) one. 
We use d models taken from Chariot and Bruzual 2007 dBruzuall 
2007a,b, hereafter CB07) to empirically define colour criteria al- 
lowing for a clean selection of z ~ 4 LBGs : 



(B-V) > 1.1 

(B-V) > 1.2 + 0.8 -(V 



H). 



In Fig.[TJwe compare the commonly adopted selection win- 
dow (left panel) to the one used here (right panel) for the CB07 
models. The blue filled squares indicate the position of a refer- 
ence LBG model with Z/Z = 0.2 and Age=100Myr from z=3 
to z=5 at 0.5 intervals, with the blue continuous lines, from left 



4 



Castellano et al.: The blue UV slopes of z~4 Lyman break galaxies 




Redshift 



Fig. 2. Upper panels: colour-colour diagrams for B-dropouts in 
the ERS (left) and HUDF (right) fields. Selected objects are 
marked in red, the blue points indicate B-dropout candidates in 
the (B-V) vs (V-Z) selection window excluded by the present 
analysis. Lower panel: combined photometric redshift probabil- 
ity distribution of our ERS B-dropout candidates (black line) 
and of galaxies selected from the (B-V) vs (V-Z) diagram (blue 
dashed line): the (V-H) selected objects have a lower probability 
of having z p hot ^ 3.5 solutions. 



to right, showing the effect of an increasing extinction from E(B- 
V)=0.0 to E(B-V)=1.0 at 0.1 intervals. Magenta and red tracks 
are for LBG models with Z/Z = 1.0 and Age=100Myr and 
300Myr respectively. Red points indicate galaxies at z=0.0-3.5 
with short star formation exponential timescales (0.1-1 Gyrs) 
and ages > 1 Gyr, while green symbols indicate constant star- 
forming models with 0.5 < E(B - V) < 1.5 in the same red- 
shift range. It appears evident that the (V-Z) colour cannot effi- 
ciently discriminate between B-dropout and lower redshift con- 
taminants, since z> 3.5 LBGs with moderate E(B-V) and both 
lower redshift dusty star-forming galaxies and passive ones oc- 
cupy the same region of the colour diagram. This effect poten- 
tially leads to a contamination in the sample of the more ex- 
tincted B-dropout galaxies, biasing the estimate of the UV slope 
distribution at z~4 and cannot be completely avoided by redefin- 
ing the relevant color criteria. The (V-H)-based criterion instead 
allows us to remove the reddest potential interlopers from the 
sample. 

Moreover, the (V-H) colour provides a cleaner LBG selec- 
tion window: a careful look at the reference LBG tracks in Fig.Q] 
shows that the (B-V)-(V-H) criteria are, in principle, unbiased 
since they are met by LBGs at z> 3.5 regardless of their E(B- 
V), age and metallicity. Specifically, there is no need to adopt a 
priori a cut in the (V-H) colour, contrary to the (V-Z) one, and 
any intrinsic constraint on the (V-H) colour due to observational 
limits corresponds to the same cut in E(B-V) for galaxies of dif- 
ferent redshifts. On the other hand, the usual (B-V)-(V-Z) cut is 
redshift-dependent, with galaxies at increasing E(B-V) meeting 



the selection c riteria at decreasing re dshifts, in agreement with 
the findings o f lBeckwifh et~aT] (120061) . 

3.2. The ERS+HUDF B-dropout sample 

Our LBG selection is naturally limited by the capability of ac- 
curately measuring large (V-H) colour terms, that is by the rela- 
tive depth of the optical and IR data-set. We limit our sample to 
S/N(H)>10cr to consider only objects with accurately measured 
colours: this choice implies that the V-band magnitude of an ob- 
ject with E(B-V)=0.5 (corresponding to (V-H)~2, i.e. /3 ~Q) is 
measured at ~3<x. On the other hand, to exclude z > 4.5 galaxies, 
we do not consider objects with S/N<3 in the V band: in prac- 
tice, given the conservative S/N cut in the H band, the latter cri- 
terion only removes one bright V-dropout galaxy from the ERS 
selection. Finally, to ensure that no AGN are present among the 
selected objects, we remove three point-like sources (Sextractor 
CLASS_STAR=0.99) from the ERS sample and check on the 4Ms 
CDFS images that none of our candidates present X-ray emis- 
sion (Fiore et al. 201 1, subm.). The final sample consists of 142 
objects at 23.9 < H < 26.5 from the ERS field, and 25 objects 
(24.8< H <27.4) from the HUDF 

We show in Fig. [2] the relevant colour-colour diagrams and 
the differences between the criteria adopted in the present pa- 
per and the usual B-dropout selection: our selected candidates 
are indicated by red points, while the blue symbols mark the 
position in the diagram of objects that would be included by a 
standard (V-Z) colour-selection, but are excluded in our sample. 
Following the previous analysis of the performance of the two 
different selections on CB07 models, the position of the latter in 
the diagrams is consistent with them being either low redshift 
interlopers or U-dropouts at z ~ 3.0-3.5. 

This can be verified through a simple test by exploiting our 
ERS photometric redshift cat alogue that also ex ploits deep K 
band and IRAC observations dSantini et al.ll2011l) . We note that 
the photometric redshifts in the ERS field have been obtained 
through the well tes ted y 1 minimization proc edure outlined in 
iGrazian et al] d2006l) and ISantini et al] d2009l) . which is based 
on a different set of synthetic spectr al templates (PEGASE 2.0 
by iFioc & Roc ca-Volmerange 1997) than those (CB07) used in 
FigQ] We determine the combined photometric redshift distribu- 
tions P(z) for the two selection criteria as the sum of the redshift 
probability distribution of the relevant B-dropout candidates. As 
shown in the lower panel of Fig|2l when objects selected through 
the (V-Z) colour are considered (blue dashed curve), P(z) dis- 
plays a lower redshift tail extending out to z ~ 2.5. On the other 
hand, the (V-H) selected sample (black continuous line) is char- 
acterized by a sharper P(z) at high redshift and has a lower prob- 
ability for secondary solutions at z < 1 : by integrating P(z) we 
estimate that the contamination in a sample selected through the 
(V-Z) colour is more than twice than in a (V-H) selected one, 
with 32% objects at z < 3.5 (14% at z < 3.0), compared to 15% 
(5% at z < 3.0) in the latter case. Photometric redshifts solutions 
at z < 1 are due to the similarity between the optical and near- 
IR SEDs of Balmer break galaxies at these redshifts, and those 
of B-dropout galaxies: in this respect, an improved colour-based 
selection is not feasible, since it would require deeper Spitzer 
data to enable a thorough study of the LBG population. 

We also exploit ERS photometric redshifts to explore dif- 
ferent criteria for both the (B-V)-(V-H) diagram and the (B-V)- 
(V-Z) one. In particular, we consider 1) 172 candidates selected 
through a less conservative (V-H) criterion, i.e. having the red- 
der selection boundary (B - V) > 1.1 + 0.7 ■ (V - H), and 2) 129 
B-dropouts isolated through a "steeper" (V-Z) colour cut aimed 
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at avoiding a redshift-dependent selection of galaxies with dif- 
ferent E(B-Vfl In the first case we find an increase of the con- 
tamination from z~ 3 - 3.5 LBGs (P(z<3.5)=18%), indicating 
that this new boundary is indeed extending towards lower val- 
ues the redshift selection window, consistently with the LBG 
tracks displayed in Fig. Q] In the second case, we still find a 
larger contamination than in our original (B-V)-(V-H) sample, 
with a 26%(12%) probability for objects at z < 3.5 (z < 3.0). 
These findings are in agreement with the tests on CB07 models 
(Sect. 13. U and show that our proposed criteria provide a more 
reliable selection for B-dropout galaxies. 

Objects having spectroscopically confirmed redshift further 
validates the efficiency of our method. Out of ten objects with 
3.5 < Zspec < 4.5 and S/N(H)>10, the (V-H) selected sample 
includes nine with 3.58 < z S pec ^ 4.4 and only fails to select a 
LBG at Zspec - 4.41. Instead, the (V-Z) colour criterion also in- 
cludes a galaxy at z sp ec = 3.11 along with those already present 
in our sample. This is also consistent with the result of intensive 
follow-up spectrocopic camp aigns of B-dropout s amples across 
all the GOODS-South field: IVanzella et alT(l2009b finds that 8 
out of 46 B-dropout galaxies initially selected through the (B- 
V)-(V-Z) diagram lie in the range 3.3 < z S n ec < 3 .5. A similar 
result is evident from Fig. 7 of IStark et a showing the 

redshift distribution of their spectroscopic sample over GOODS- 
North. We show in the next section and in Fig.|6]how the inclu- 
sion of these potential interlopers in our sample would affect the 
determination of the UV slope distribution. 



4. Estimating UV slopes and rest-frame magnitudes 

We adopt the common power-law approximation for the UV 
spectral range oc and we estimate the slope /3 by fitting 
a linear relation through the observed I, Z, Y, J AB magnitudes 
of the objects: 



Mi = -2.5 • 06 + 2.0) • log(Ad + c 



(1) 



where M, is the magnitude in the i-th filter at effective wave- 
length Aj. The filters we use span the rest-frame wavelength 
range A =* 1500 - 2500 A at z ~ 4, providin g a /? estimate whic h 
is consistent with the original definition by Calze tti et al.l dl994h . 
and that is independent on the exact redshift of the sources, 
which only affects the intercept c in Eq. Q] We underline that 
such estimate of ft is unbiased with respect to the colour selec- 
tion being based on a different set of filters. We present in Fig. [3] 
four examples taken from the ERS sample, along with a compar- 
ison between the UV slopes determined by this method and the 
relevant best-fit SEDs. 

The large wavelength baseline we adopt gives ft estimates 
that are significantly more stable than e stimates based on the I—Z 
colour alone (e.g. lBouwens et aill2009l) : as an example, we show 
in Fig. [4] a comparison of the ft uncertainty distribution of the 
two methods for H~27 objects from the HUDF mock catalogues 
presented in Sect. 15.11 

Fig. [5] shows the resulting UV slope ft for the objects in our 
sample as a function of their H band magnitude: most of the 
galaxies have blue UV slopes (ft ~ -2), apart from a small 
number of red LBGs mainly populating the brightest end of the 
diagram. We explore how the distribution of ft would change 
when adding to our sample objects that are included by a stan- 
dard B-dropout selection diagram, but that are not selected when 
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Fig. 3. A comparison of the fitted UV slope (red continuous line) 
and best-fit SED of four objects in our sample having best-fit 
photometric redshift in the range z P hot = 3.5 - 4. 1 . The observed 
magnitudes (filled squares and errorbars) are, from left to right: 
B, V, I, Z, Y, J, H. The linear-fit is performed on I, Z, Y, J. The 
lower and upper horizontal axis show the observer-frame and the 
rest-frame wavelength, respectively. 
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1 defined on the basis of CB07 LBG tracks: (B-V) > 1.1 A (B-V) > 
1.1 + 1.5-(V-Z)A(V-Z)< 1.6 



Fig. 4. Comparison of uncertainties on ft derived from a linear 
fit to the observed I, Z, Y, J magnitudes (histogram) and from 
the (I-Z) colour only (red curve) for simulated objects with ob- 
served UV slope j6/„. The real UV slope ft mput of each model is 
measured on its unperturbed magnitudes. 



their (V-H) colour is considered. We perform this test on the 
ERS to explore also their photometric redshift distribution. As 
shown in Fig. [6] these candidates (magenta dashed histogram) 
are on average redder than those already present in our sam- 
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Fig. 5. The UV slope /3 of the ERS (open triangles) and HUDF 
(filled squares) B-dropout objects as a function of their observed 
H band magnitude. Errorbars are given by the uncertainty on the 
linear fit on their observed I, Z, Y, J AB magnitudes. 
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Fig. 6. Distribution of the UV slope values from the ERS sample 
(black hollow histogram). The magenta dashed histogram rep- 
resent the objects that are not included in our sample but are 
selected through a standard (V-Z) colour-cut: most of them have 
Zphot < 3.5 (blue filled histogram). The blue arrow marks the 
position of a z spe c = 3.11 object included in this last sample. 

pie, with 38 out of a total of 47 having a best-fit photometric 
redshift z p hot < 3.5 (blue filled histogram). Moreover, the only 
spectroscopically confirmed object at z spec = 3.11 that we ex- 
clude through the (V-H) criterion has a very red UV slope with 
P - -0.75 (blue arrow in Fig. [6}- Notably, all the reddest objects 
(/J ~ 0) of the (V-Z) selected sample can be considered inter- 
lopers on the basis of their photo-z. These results are consistent 
with the analysis of LBG models presented in Sect. 13. II showing 
that the usual (B-V)-(V-Z) cut preferentially selects objects with 
large E(B-V) at z -3-3.5. 

A more detailed analysis of the LBG population can be given 
by looking at how UV slopes change as a function of the rest- 
frame UV magnitude: to this aim, a linear fit of the UV SED 
gives us a further advantage, namely we can interpolate the slope 
of each object to estimate the rest-frame magnitude at 1600A. 
To provide an estimate for Mi6oo> we consider all objects as lo- 
cated at the same redshift, which for simplicity we assume to be 
z — 4: the effects of this assumption are self-consistenly taken 
into account in the statistical analysis discussed in the following 
sections. 

5. A statistical analysis of the UV slope distribution 

5.1. Uncertainties and systematics in p and M mm estimates 

We first start with a set of Monte Carlo simulations reproducing 
the overall properties of our observations. We use the CB07 li- 
brary to produce a large set of simulated galaxies at 3.0 < z < 5.0 



with expected magnitudes in the same filter set as in the ERS 
and HUDF observations. We consider models with constant star- 
formation histories and within the following range of physical 
parameters: Z/Z = 0.02, 0.2, 1 .0, 0.01 < E(B - V) < 1.0, 
0.01 < Age < lGy r. We assume a lSalpeterl d 19551) IMF and a 
ICalzetti et alj (2000) attenuation law, while the transmission of 
the IGM is treated according to lFan et all (12006). We exploit this 
library to generate mock catalogues of 5 x 10 5 objects for each 
field by perturbing magnitudes in all the bands to match the rel- 
evant depths in the observed ones. Our simulations are designed 
to reproduce both the average and the scatter of the S/N distri- 
bution as a function of magnitude to provide the closest match 
between observed and simulated data. The simulated catalogues 
are then analysed just as the real ones to determine the obser- 
vational effects on galaxies of different magnitudes and spectral 
slopes. 

Our Monte Carlo approach, while being much less time con- 
suming than two-dimensional imaging simulations, allows us to 
estimate the effect of magnitude and colour uncertainties in the 
dropout selection and UV slope analysis. A further advantage of 
such a purely empirical approach is that it is free of assumptions 
on the typical size and size-luminosit y relation of the LBGs, 

which, as shown in iGrazian et al. i (HoHb 

are a critical factor to 
estimate the completeness correction for faint objects. We any- 
way randomly select a subsample of 20000 model LBGs from 
the same CB07 sample and perform a simple two-dimensional 
simulation to compare its results to those obtained through our 
empircal method. We adopt as morphology template the stack- 
ing of ten bright candidates from the ERS field, and we perform 
various runs inserting 200 objects in the images each time af- 
ter masking the regions where real objects were detected. The 
inserted objects are then recovered and analysed with the same 
procedure used for the real catalogues. A comparison between 
simulated objects obtained in this way and the mock catalogues 
discussed above shows a good consistency, therefore we will 
base our analysis on the latter ones that allow us to perform 
larger and less time-consuming simulations. 

In the uppermost panel of Fig. [7] we show for both fields 
the detection completeness at different observed fluxes of galax- 
ies with a given spectral slope, as indicated by the fi mpul of the 
model measured on its real unperturbed magnitudes. For the 
ERS sample the completeness curves remain approximately flat 
for Pinpu, < -1.5 and rapidly decrease at fi mpul > -1, especially 
for faint objects. The completeness in the HUDF is higher be- 
cause of the deeper optical images: it remains close to one up to 
flinput 55 -0.5 for bright objects, while fainter ones are still effi- 
ciently selected as long as /3i„ pu , < -1.0. The combination of the 
two fields thus allows us to sample a wide range of UV slopes 
and magnitudes. 

In Fig. |71 we also compare input and measured values of 
M1600 and [3 for galaxies taken from our simulations. In partic- 
ular, we explore the dependence of uncertainties and systemat- 
ics on the spectral shape of the galaxy, as indicated by /?;„ piJ f, 
and on the observed magnitude. The latter test is of particular 
importance since we are analysing the colour distribution of a 
flux-limited sample. The uncertainties on both M1600 and /3 are 
relatively small, with cr M1 6oo — 0.1 - 0.2 and crp - 0.1 - 0.3. 
However, they become larger for faint and red objects in the 
sample with respect to blue/luminous ones. The measured UV 
slope Pji, does not show any large systematic offset as long as 
Pinpui < -1, a small trend towards underestimated values ap- 
pears for the faintest red objects in the ERS sample. This is due 
to the shallower optical data in that field. The estimated UV 
rest-frame magnitude M1600 shows a negligible offset of A ~ 
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Fig. 7. From top to bottom: detection completeness and differ- 
ence between output and input j3 and Mi6oo from our ERS (left) 
and HUDF (right) mock catalogues. All quantities are computed 
for objects in three observed magnitude bins centred at: H=24.5 
(ERS) and H=25.5 (HUDF) (black continuous lines and filled 
squares); H=25.25 (ERS) and H=26.3 (HUDF) (dashed blue 
lines and open squares); H=26.0 (ERS) and H=27.0 (HUDF) 
(red dot-dashed lines and crosses). 



0.05-0. 1 mags with respect to the real value for f3i nput < - 1 with 
redder objects being underestimated by A ~ 0. 15-0.2 mags. The 
conservative selection S/N(H)> 10 helps keeping the uncertain- 
ties small: a similar test for objects of S/N(H)~ 5-10 shows 
that the uncertainty on the UV slope is as large as 0.6 for faint 
galaxies. 



5.2. Reconstructing the UV slope distribution: the method 

The tests discussed in Sect. 15. II show that the completeness and 
the accuracy on both /3 and UV rest-frame magnitude estimates 
depend on the observed magnitude as well as on the colour of 
the galaxies. This implies that any measured UV slope distribu- 
tions at different luminosites arise from the underlying intrinsic 
one after being convolved by a combination of observational ef- 
fects. To estimate the dependence of /? on UV magnitude it is 
therefore not sufficient to apply simple corrections on the ob- 
served distributions but it is necessary to adopt a more com- 
prehensive approach. Such an approach has to be based on as- 
sumptions on the properties of the intrinsic UV slope distribu- 
tion, likewise any "de-convolution" procedure. We apply here a 
maximum likelihood technique which is analogous to the one 
used to estimate the galaxy Luminosity Function (LF) where a 
Schechter LF is assumed and its best-fit parameters are obtained 
by comparing simulated number counts to the observed ones. We 
start from the assumption that the probability distribution func- 
tion of UV slopes in a given magnitude range, PDF(J3), follows 



a given functional form, whose average and standard deviation 
are to be determined. For any given set of parameters, and for 
each magnitude bin, we extract from our CB07 libraries 50000 
objects reproducing the given distribution and perturb them with 
noise as discussed in Sect. 15.11 to introduce the effect of ob- 
servational uncertainties and biases. Finally, we compare the 
output distribution of /? values to the observed one in each bin 
through a maximum likeliho od estimator £, (e.g. lBouwens et all 
I2008t ' ICastellano et al.ll20lob and minimize A^ 2 = -2.0 • ln(£). 
Unfortunately, while a Schechter form in the case of the lumi- 
nosity function is suggested by theoretical arguments and has 
been verified through many samples at low and high redshift, 
the present work constitutes the first attempt of a parametri- 
cal reconstruction of the UV slope distribution PDF(/?). We are 
thus forced to start from reasonable hypothesis on the analytic 
shape of PDF(y6), as suggested by our data and by previous non- 
parametric analysis. We adopt two different analytic PDF(y6) as 
working hypothesis. We first assume a Gaussian distribution for 
the [5 values. We note that such a distributi on has already bee n 
suggested by other authors in the past (e.g. lBouwens et aiT2 009). 
As an alternative, we assume a log-normal distribution PDF(yS) 
with (3 = -2.5 as lowest, asymptotic value. This is motivated 
by the asymmetric distribution observed in the brightest mag- 
nitude bins (Fig. [8} and corresponds to the hypothesis that the 
high-z galaxy population is made of a bulk of low-extincted star- 
forming LBGs along with a number of rare objects at increasing 
E(B-V) values. 

We change the average and dispersion of both input distri- 
butions across a 25x25 grid in the range (J3) = -2.5, -1.5 and 
r.m.s.(/?)=0.03,0.7. 

5.3. Results 

In order to analyse the relation between magnitude and UV 
slope, we divide the sample into four M1600 bins (Fig. [8). This 
choice is somewhat arbitrary and is motivated by the compet- 
ing needs of having magnitude bins as narrow as possible, but 
with a sufficient number of objects allowing the statistical anal- 
ysis described in the following. In practice, given the different 
magnitude ranges probed by the two fields, the ERS objects will 
be used to constrain the first and second bin, the third one will 
be analysed by combining data from both fields, while HUDF 
objects will be used for the faintest magnitude bin. 

The resulting best-fit distributions (Table []]) indicate that 
most of the B-dropouts have blue UV slopes: in particular, at L < 
L* galaxies tend to be very blue regardless of their M1600, while 
the brightest objects only (first bin) are, on average, slightly red- 
der than faint ones. Indeed, the best-fit Gaussians are remarkably 
similar in the three magnitude bins at M1600 > -21 .0 with an av- 
erage of (J3) -2.1 and an intrinsic r.m.s. ^ 0.3, while brighter 
galaxies have a best-fit Gaussian with (J3) - -1 .9, but similar 
dispersion. Given the low uncertainties, this small difference be- 
tween bright and faint galaxies is significant at ~ 2<x. The log- 
normal fit yields similar results, with bright galaxies (first bin) 
redder than Migoo > -21.0 objects, the only difference being 
that in this case, given the asymmetric shape of the distribution, 

the bulk of the objects would have /3 2.0 in the first bin and 

/3 ~ -2.2 at M1600 > -21.0. We verified that these results are 
robust against small changes in the colour selection criteria: ob- 
jects selected with the redder boundary in the (B-V) vs. (V-H) 
diagram discussed in Sect. 13.21 yields best-fit values within the 
68% c.l. region given by our "fiducial" sample. 

The binned data and best-fit distributions along with the 
maximum likelihood contours for the parameters are shown in 
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Fig. 8. Maximum likelihood contours (68%, 95% and 99% c.l. in green, blue, and red respectively) for the Gaussian and lognormal fit 
at different magnitudes. The histograms are the observed data. The red continuous line in each bin is the relevant best-fit distribution. 
The blue line show how it appears after being "convolved" by observational effects. 



Fig. [8] the need for an approach like the one carried out here can 
be appreciated by looking at how the input best-fit distributions 
in the various magnitude bins (red curves in Fig. [8) are modified 
when they are self-consistently "convolved" with observational 
effects through simulations (blue curves). 

Finally, Fig. [9] shows the distribution of observed points in 
the /3 - M1600 diagram along with the average and intrinsic dis- 
persion of the best-fit Gaussian and lognormal distribution. 

While our tests consistently indicate that most of the B- 
dropouts have blue UV slopes, it does not allow us to determine 
which is the functional form that best reproduces the data. To 
this aim, we compared the observed UV slope distributions to the 



best-fit Gaussian and log-normal PDF(/3) in each bin through a 
Kolmogorov-Smirnov test. We find a significant difference only 
in the first bin where the Gaussian fit yields a low KS probabil- 
ity (<50%) because of the presence of a tail with few red objects 
at /3 ^ -1.0. In the remaining bins the KS test does not allow 
us to discriminate between the two distributions. It is interesting 
to note that a constant distribution of /3 values within each bin 
is completely excluded by our maximum likelihood test at all 
magnitudes. However, we cannot exclude that the distribution is 
bimodal, with a "red" peak at values of /3 larger than accessible 
to Lyman Break selections, or that it is a more complicated one, 
as it would be the case of Gaussians with non-standard kurto- 
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Fig. 9. The ft - Mieoo diagram for our B-dropout sample. Open triangles and filled squares represent individual objects in the ERS 
and HUDF fields, respectively. Errorbars on the observed points are given by the uncertainty on the linear fit on their observed 
magnitudes. Blue and red circles and lines indicate the result of the statistical analysis described in Sect. [5] which takes into account 
systematic effects. Blue open circles and errorbars indicate the average (J3) and intrinsic r.m.s of the best-fit Gaussian distributions 
in each of the four bins: uncertainties on (fi) are much lower and they are not visible in this scale. Red filled circles for the best-fit 
log-normal distributions. The Gaussians best-fit values are plotted at the average magnitude of the observed points in the relevant 
bins, log-normal ones with a small magnitude offset for clarity. The green continuous line is the /3 — UV relation by Bou wens et al.l 
(2009). We show for comparison a linear fit on the Gaussian and log-normal best-fit points (blue dashed and red continous lines): 
we underline however that our analysis indicates a significant difference only between the first and second magnitude bins. 



Table 1. UV slope distributions: best-fit parameters and 68% c.l. 
uncertainties. 
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sis and/or skewness. While the first case can only be probed by 
combining LBG sele cted samples with far -IR and submillime- 
tre selected ones (e.g. lWardlow et a0l201 ll) . the latter would re- 
quire fitting more parameters than it is meaningful to do with the 
presently available data. 

The physical interpretation of these results is not straightfor- 
ward, since several factors determine the UV slope: dust redden- 
ing, varying ages, metallicity, IMF, and star f ormation histories 
(e.g. lBouwens et al.ll2009HWilkins et al.ll201 lh . However, the ef- 
fect of dust is by far the most relevant and with the presently 
available data it is not possible to constrain how much other 
properties contribute to the variation of the UV slope distribution 
with magnitude, or if they only affect the intrinsic dispersion of 
the PDF(y6) distributions. For this reason, and to provide com- 
parison with other works, in the following we will interpret our 
results by directly converting the UV slope into dust extinction. 

In this context, our results imply that the average extinction 
at 1600A (following iMeurer et alJfl999h is A 160 o - 0.65 in the 
firs t bin and A1600 - .25 in the faintest bins. Deviations from 
the Meurer et al. (1999) relation have been found in z>4 QSOs 
by iGallerani et al. (2010): their MEC attenuation curve would 
imply a higher attenuation, A1600 - 1-4 and A1600 - 0.5 for 
bright and faint galaxies respectively. 



Previous determinations of the UV ex ti nction of z ~ 
have be en provided b y lOuchi et al. d2004l). 
d2006l). [Overzier et alJ d2008l) iBouwensetal 
Lee etalJ (120111). 



4 LBGs 
jeck with et al.l 



(2009) and 



lOuchi et all (120041) found an average E(B-V)=0.15 + 0.03 
for their M < M* B-dropouts and no dependence on the ap- 
parent (i.e. uncorrected) luminosity within the magnitude range 
they probe. Our estimate at similar magnitudes (Mi^oo < -21.0) 
is lower: when converting the relevant best-fit Gaussian PDF(/3) 
into a corresponding E(B-V ) distribution following the equa- 
tions in Calze tti et al.l d2000l) . we obtain an ave rage (E(B - V)} = 
0.05^. On the other hand. iLee etaP (1201 lh find a strong de- 
pendence of fi on UV magnitude in a large sample of M < M* 
z ~ 3.7 galaxies from the NOAO Deep Wide Field Survey. 
Although they probe much brighter magnitudes than those avail- 
able in our sample, a comparison can be made between their 
faintest bin centred at Mnoo = -21.43 and the brightest galax- 
ies from our ERS sample. They find /3 = -1.78 + 0.28 at these 
magnitudes which falls within the 68% c.l. region given by our 
ma ximum li kelihoo d test. 

iBeckwith et alJ d2006l) find blue slopes for both bright and 
faint z ~4 galaxies in the HUDF, suggesting that they are affected 
by little du st extinction. A quantit ative comparison is not possi- 
ble because Be ckwith et al.l d2006) do not explicitly compute the 
value of the UV slope in the different samples. However, their 
result is in agreement with our findings of blue UV slopes in all 
L < L„ magnitude bins, w hich is the range probed by the HUDF. 

lOverzier et al.l d2008l) provides an estimate of the UV slope 
of LBGs and Lyman Alpha Emitters (LAE) in a z~4. 1 proto- 
cluster finding a global average (/3)=-1.95 which is in agreement 
with our results. 

On the other hand, our results are markedly different from 
those obtained bv lBouwens et alJ d2009l) . They find redder j3 av- 
erages and a steep correlation between f3 and UV magnitude up 
to very faint luminosities: their relation is shown as a green line 
in Fig. [9] for comparison we show as blue (red) lines a simple 
linear regression of our best-fit Gaussian (log-normal) average 
slope (fi) as a function of Mi6oo- 
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These discrepancies can in part be explained by the differ- 
ences between various approaches, especially by the different 
colour selections as discussed in Sect [3] Specifically, we ex- 
plo red whether the g reates t differences between our analysis and 
the iBouwens et al.l d2009l) one can be responsible of the dis- 
crepancies: we mimicked their approach by adopting the B-V, 
V-Z criterion (i.e. including in our sample the objects shown in 
Fig. O and we adopted the same definitions for/? and Muv based 
on the I and Z observed magnitudes. With the inclusion of the ad- 
ditional red candidates, at bright magnitudes we obtain average 
UV slopes redder (J3 ~-1.7) than in the case when potential in- 
terlopers are excluded through their (V-H) colour. However even 
wi th this inclusion, the v alues still are not as red as those found 
bv lBouwens et al.l d2009lPl 

A completely different approach in the estimate of dust ex- 
tinction in high redshift galaxies is provided by stacking analy- 
sis at radio wavelengths to determine their average un corrected 
SFR. I n this respect, the analysis of z ~ 4 LBGs by iHo etall 
based on VLA 1 .4GHz observations, favors lower val- 
ues for dust extinction than those inferred from past UV analy- 
sis, in agreement with our findings. However, a comparison be- 
tween the two approaches is c omplicated by seve ral factors and 
assumptions (see discussion in lCariffi et alJ l2008). 



6. The dust corrected Star Formation Rate Density 

The parametric estimate of the UV slope distribution at different 
magnitudes gives us the remarkable opportunity of computing 
the SFRD in a straightforward way through the following equa- 
tion: 

SFRD = - l '^ 27 j dL j dA- PDF {A, L) ■ 10° 4 A • L ■ <D(L)(2) 

where the constant factor is from Mad au et al.l (Il998h . 
A=Ai6oo sucn that PD F(A,L) is un iv ocally related to the 
PDF(J3,L) through the IMeurer et ail ( 119991) equation, and 
<£(L) is the UV luminosity function at 1600A. For the latter 
we ad opt the best-fit Schechter parameters by IBouwens et alj 
(2007|) for consistency with previous results. In practice, we 
compute the integral in each luminosity bin by consider- 
ing the relevant best-fit distributions PDF(/3), and extrapolat- 
ing the last one up to M1600 = -17.5 (corresponding to 
L = 0.04L*). When the Gaussian PDFQ3) are considered 
we obtain log(SFRD)=-1.093+° j^ Me/yr/Mpc 3 , taking into 
account both the 68% c.l. uncertainty on each PDF(/3) and 
the uncertainty on the LF parameters. Using the best-fit log- 
normal PDF(j0) in Eq. [2] we obtain a slightly higher value 
log(SFRD)=- 1.086+°^ M Q /yr/Mpc } because of the larger ef- 
fective correction for dust e xtincion. Our SFRD e stimates are 
~55% lower than the one by IBouwens et alj d2009l) : this not a 
surprise since we find bluer /3 at bright and intermediate magni- 
tudes, hence we apply a lower correction for dust. Interestingly, 
our SFRD values are c onsistent with the estimate obtained by 
I Thom pson et al. (2006) on a photometric-redshift selected z ~4 
sample and through the use of SED-fitting to estimate extinction 
correction. This is in agreement with the discussion in Sect. [3] 
and [4] on the differences between our H-based selection and the 
usual optical-based one. 

1 We note that Bouw ens et alj d2011l) have recently presented a re- 
fined estimate of the average UV slopes where WFC3 IR information 
is also exploited, as in the present paper. Their findings show a good 
agreement with our estimates. 



The largest uncertainties in the SFRD estimate are given by 
our lack of knowledge of the attenuation curve at these redshifts, 
and by the possible presence of highly-obscured star-forming 
objects, like the sub millimetre galaxies (S MGs) well studied at 
lower redshifts (e.g. IChapman et al.ll2005|). As a n example, the 
MEC attenuation curve bv lGallerani et alj (2010), which is flat- 
ter than the Calzetti et al. one in the UV, yields a SFRD nearly 
three times larger: log(SFRD)^ -0.65 M Q /yr/Mpc 3 . On the 
other hand, a significant contribution from z ~ 4 SMGs (e.g. 
Daddi et al. 2009) would imply higher values of the SFRD than 
those measured from LBGs alone. 



7. Dust extinction as a function of intrinsic UV 
luminosity 

Galaxy evolutionary models suggest that UV bright galaxies 
tend to be dustier than faint ones, as a natural outcome of the re- 
lation between SFR, mass and metallicity. Our findings of a red- 
der UV slope distribution at bright magnitudes apparently con- 
firm this prediction. However, such analysis shows how galaxies 
of different UV slopes (which is a proxy of their dust content) are 
distributed as a function of their extincted UV magnitude. The 
effect of extinction is thus present on both quantities and a clear 
characterization of how the dust content varies across the galaxy 
population is not possible: the usual j6-UV magnitude relation, 
as presented in Sect. [5] is mainly useful to estimate a correc- 
tion to the observed SFRD, as discussed above, but do not offer 
a straightforward insight on the physical properties of Lyman 
break galaxies. 

To provide a clearer view on this issue we show in Fig. [10] 
the relation between f3 and the dust-corrected UV magnitude 
Midoocorr = -^1600 - Ai6o(). We also show the conversion of 
both in more grou nded phys ical quantities: Aim o (following 
IMeurer et all 1 19991) and SFR dMadau et al .11 19981) respectively. 
Naturally, degeneracies between extinction and other physical 
parameters can bias this direct conversion, and quantities are not 
completely independent, since the extinction correction is itself 
related to the (5 of the single objects. Nonetheless, the simplic- 
ity of our analysis provides a more straightforward view on such 
observational degeneracies than in the case of common SED- 
fitting estimates. We also plot in Fig. [10] the 30% and 90% de- 
tection completeness levels, as estimated from our simulations, 
for galaxies of the two fields. 

From this plot it appears evident that at large intrinsic lumi- 
nosities (Mi600co/r 55 -23.0, corresponding to SFR> 80 M Q /yr), 
only dusty star-forming galaxies are found. Once the incom- 
pleteness of our selection is taken into account we obtain for 
objects with SFR> 80 M G /yr a surface density of 0.12 arcmirT 2 
when 2 < Ai6oo < 3 and an upper limit of 0.02 arcmirT 1 for less 
extincted galaxies. 

Among galaxies with lower SFR, we detect many with a 
much lower amount of reddening. As a result a trend of decreas- 
ing extinction as a function of corrected UV luminosity seems to 
be in place. However, this effect may be entirely a systematic ef- 
fect due to the limits of our observations (as shown by the lines 
in Fig. [Tol l that do not allow us to ascertain whether this trend 
is real or it is due to selection effects even at moderately high 
luminosities. 



8. Summary and Conclusions 

In this paper we have presented a new approach to the selection 
and to the analysis of the UV slope distribution of z~4 Lyman 
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Fig. 10. The relation between (3 and the dust-corrected UV magnitude for ERS (open triangles) and HUDF objec ts (filled squares) in 
our sampl e. The upper and th e right axes show the conversion of unextincted magnitude and UV slope in SFR (Ma dau et al.ll 19981) 
and A i6oo (M eurer et al.ll 19991) respectively. The red continuous lines and the blue dashed ones indicate the 90% and 30% selection 
completeness in the two fields. 



Break galaxies. Thanks to the deep WFC3 near infrared obser- 
vations of the ERS and HUDF fields, combined with the existing 
HST-ACS optical data, we perform a B-dropout selection that is 
based on a (B-V) vs. (V-H) diagram on H-detected samples of 
both fields. Following the analysis of galaxy spectral libraries 
and of our own photometric redshift catalogue on the ERS, we 
have shown that this selection allows us to efficiently exclude 
low redshift interlopers and dusty z>3 galaxies that might be 
present in B-dropout samples selected through optical colours. 
The sample we select on the basis of the newly defined colour 
criteria comprises 142 and 25 objects at S/N(H)>10 in the ERS 
and HUDF fields respectively. We estimate the UV slope of our 
B-dropout candidates through a linear fit on their I, Z, Y, J mag- 
nitudes, and we interpolate the same fit to provide an estimate 
of their UV rest-frame magnitude after assuming that all sources 
are located at z=4. 

Despite the conservative selection and the greater accuracy 
in the estimate of /3 provided by a linear fit, we show that uncer- 
tainties and systematics arise when measuring the colour distri- 
bution of such flux-limited samples. However, the effect of ob- 
servational uncertainties can be effectively taken into account 
by adopting a parametric approach combined with detailed sim- 
ulations, and the estimate of the typical UV slope at different 
magnitudes can be performed following hypothesis on the shape 
of the undelying intrinsic distribution. We carry out a parametric 
maximum likelihood analysis by assuming that the distribution 
of ft is either Gaussian or lognormal. We allow the relevant pa- 
rameters of the two distributions to vary across a wide range of 
values and find the best-fit values by comparing the simulated 
and observed counts through a maximum-likelihood estimator. 

We find that z ~4 LBGs are characterized by blue UV slopes, 
suggesting a low dust extinction. The best-fit Gaussians are re- 
markably similar in all the magnitude bins at Mi6oo > -21.0 
with an average of {ft) =t -2.1 and an intrinsic r.m.s. =t 0.3, 
while brighter galaxies have a best-fit Gaussian with slightly 



redder average {(ft) =* -1.9) but similar dispersion. A similar 
result is obtained when log-normal distributions are considered, 
the main difference being that in this case, given the asymmet- 
ric shape of the distribution, the bulk of the objects would have 

ft ~ -2.0 in the first bin and ft 2.2 at M 16 oo Z -21.0. 

When interepreted in term of varying dust extinction, our results 
imply that z~4 LBGs have a low dust content, with a typical 
A1600 - 0.65 at L < L * and A 1600 - 0.25 for faintest galaxies, 
following lMeurer et all dl999l) . 

Ou r results are in general agr eement with the prev ious find- 
ings of lBeckwith et alJ (l2006l) and lOverzier et al . (2008), and are 
consistent within the uncert ainties with the result on L ~ L* 
galaxies bv lLee et alJ d201 ll). On the other han d our results dif- 
fer from thos e obtained by Ouchi et all d2004l) for bright galax- 
ies, and from lBouwens et all d2009l) that find redder ft averages 
and a steeper correlation between ft and UV magnitude. We find 
that these discrepancies are only in part due to the differences 
in the selection (IR-based vs. optically-based) and in the analy- 
sis of the B-dropout samples, since such red averages for the UV 
slope are not found in our sample even in the case when potential 
interlopers selected through their optical colours are included. 

The low dust extinction we find points to a SFRD in the 
LBG population that is significantly revised downward with re- 
spect to previous estimates. If we compute the corrected SFRD 
from the UV LF by exploiting our parametric estimate (Eq. 13, 
we obtain log( S FRD) ~ -1.09, which is ~55% lower than the 
Bouwens et al ] d2009h estimate. However, this result is based 
on the assumption of a C alzetti et al. attenuation law: when the 
MEC attenuation curve of iGallerani et al.l d2010h is assumed we 
obtain a nearly three times higher SFRD at z ~4. 

Finally, we discuss how the UV slope of z~4 galaxies 
changes as a function of the dust-corrected UV magnitude, i.e. 
as a function of the SFR. Because of the difficulties in se- 
lecting highly extincted objects, we show that constraints can 
only be given on the number density of UV bright galaxies 
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(SFR> 80 M B lyf). We find that highly extincted galaxies with 
2 < A i(,oo < 3 are much more common (surface density of 
0.12 arcmirT 2 ) than dust-free ones (< 0.02 arcmirT 2 ). This im- 
plies that most galaxies with a high SFR (> 80 M /yr) are highly 
extincted objects. Among galaxies with lower SFR, we detect 
many with a much lower amount of reddening, although current 
observational limits prevent us from detecting those with high 
exctinction, if they exist. 

Although the now available WFC3 observations allow for a 
refined selection of z~4 samples, improvements are needed to 
further constrain the dust properties at high redshifts. Deeper 
spectroscopic follow-up observations are necessary to validate 
both IR and pure optical colour selections. Such samples will 
also give the opportunity of complementing the study of UV 
slope distributions with the analysis of the Lyo' EW distribu- 
tion whic h might also be relat ed to the dust content of galax- 
ies (e.g. iPentericci et alj 120071) . Larger near IR data-samples, 
as those that will be available in the near future thanks to the 
CANDELS survey, will allow us to deepen this kind of analysis. 
While CANDELS/Wide fields will enable a comparison of the 
UV slope distribution of IR-selected LBGs to the past optical- 
based analysis of large L>> L* samples, the CANDELS/Deep 
observations will allow us to strengthen statistical constraints at 
faint/intermediate magnitudes. 

Instead, a thorough characterization of the high-redshift pop- 
ulation and a comprehensive analysis of the relation between 
SFR and dust content will require future instruments both in the 
optical and in the IR capable of observing highly extincted ob- 
jects with SFR of few solar masses per year. 
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